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Abstract 



We study some properties of the deformed Sine Gordon models. These models, presented by 
[■"-« I Bazeia et al, are natural generalisations of the Sine Gordon models in (1+1) dimensions. There 

are two classes of them, each dependent on a parameter n. For special values of this parameter the 
models reduce to the Sine Gordon one; for other values of n they can be considered as generalisations 
^^ I of this model. The models are topological and possess one kink solutions. Here we investigate the 

" ' existence of other solutions of these models - such as breathers. The work is numerical and we find 

that the breathers, as such, probably do not exist. However, we show that some of these models, 
namely, the n = 1 of the first class possess breather-like solutions which are quasi-stable; ie these 
^..^ ^ 'quasi-breathers' exist for long periods of time (thousands of periods of oscillations). These results 

l/^ ' are found to be independent of the discretisation used in the numerical part of our work. 



PACS numbers: 05.45.Ac, 05.60.Cd 



oo ! 1 Introduction 

O 

Recently [1] Bazeia et al have presented two classes of models in (1+1) dimensions which are natural 
r> I generalisations of the well known Sine Gordon model. The models in each class depend on one 

parameter n. For special values of this parameter, namely n = 2 for the first class and n = 1 for the 
second class, the models reduce to the Sine Gordon one; for other values of n they can be considered 
as generalisations of this model. The models are defined by the Lagrangian 

'^ " (1) 
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When n = 2 for the first class, and n = 1 for the second, the models reduce to the Sine Gordon 
model (for ip = (p ox 2(j) respectively). 

The models are topological and for any n they posses one kink solutions. 
For the first class of models these one kink solutions are given by 

exp (27 ( X — xo — ut)) 



(p(x, t) = sin ^ . , , , 

^^ ' ' Ll + exp(27(x-xo- nt))_ 

These kink solutions satisfy the usual Sine- Gordon boundary conditions, namely: 



(4) 



TT 

(/9 ( X ^ —00) — > 0, ip{x ^> 00) > — (5) 

The kink solutions of the second class are given by 

ip{x^t)= 2tan^^ {exp^-j {x — xo — ut))) ". (6) 

Their boundary conditions are 

if{x^ —00) > 0, (/? ( X ^ 00) > TT" (7) 



In these expressions 7 is the usual relativistic factor , ^ , where u is the velocity of the kink. 

The kink solutions have been studied in the original paper of Bazeia et al [1] and in a recent paper 
by al Alawi et al [2]. These papers showed that as these models are topological the kink solutions are 
stable. When the kink solutions were sent towards various obstructions, in the form of a barrier or 
a potential hole, the behaviour of the kinks resembled the behaviour of the kinks of the Sine Gordon 
model; ie the kinks could get trapped in the holes and at some very precise velocities could even get 
reflected by the holes. 

When one considered the dependence on n of these properties (like the critical velocity below which 
the kinks were trapped) it was clear that the Sine Gordon model was not special from this point of 
view; for example, the critical velocity for the transmission was lowest for n = 3, and not n = 2, in 
the first class of models. 

Hence it is natural to check whether the models of Bazeia et al share other properties of the Sine 
Gordon model, like integrability and/or the existence of other fine energy solutions. Of these finite 
energy solutions the easiest to study are the multi-kink solutions and breathers. 

Unfortunately, so far we have failed to find such solutions analytically. So we have been forced to 
look at them numerically. 

In this paper we report some results of our numerical studies. In the next sections we present 
some results of our studies of breathers. These results were obtained in the usual discretisation. Our 
simulations were performed using 4th order Runge-Kutta method of simulating time evolution. The 
spatial lattice streched from —50 to 50 and the lattice spacing was dx = 0.01. The time evolution 
spacing dt was taken as dt = 0.0001. 

Later we look at the topological discretisation of Speight and Ward [3] and apply it to our case. We 
find that that the results are not very different from what we have seen using the usual discretisation 
so that the fact that our 'breather-like' states are long-lived but not really stable is not a numerical 
artefact. In the following section we say a few words about two kink states. The final section contains 
some further comments and a discussion of our future plans. 



2 Breathers 

In the Sine Gordon model (ie for n = 2) the potential is given by 

Vi^) = ^sm^lip) (8) 

and one has analytical expressions for breathers. For n ^ 2 the equations of motions are more 
complicated and, so far, we have not found analytically any breathers. However, as is well known 
[4], breathers are bound states of kinks and antikinks and, as such, can be also found by taking a 
configuration consisting of a kink and antikink (not too far and not too close to each other) and 
letting it evolve with time. The field quickly readjusts itself to that of a breather with some small 
deformations which move out towards the boundaries. These deformations represent an extra energy 
which is being sent towards the boundaries. We absorb it and the field becomes that of the breather. 

In fig. 1 we present the energy of three such evolutions for the case of a kink/antikink configuration 
in the Sine Gordon model. We generate such a configuration by taking our field (p corresponding to 
a kink located at x = —d for x < and the corresponding antikink (located at x = d) for x > 0. 
For d > 4.5 the resultant configuration is quite smooth around x = 0, but to make it even smoother 
we have taken the average values of kink and antikink for —0.5 < x < 0.5. We have checked that 
the results do not depend on this 'smoothing' procedure. Here we present the results of three such 
simulations. In fig la - the present the total energy for the case of d = 3, in fig lb - d = 5 and in fig 
Ic of d = 8. We note only a small initial decrease in energy after which the energy is well conserved. 
This initial decrease is associated with the rearrangement of the field (from our initial configuration 
to that of the breather one) and then the field evolves as a genuine breather. In fig 2 we present 
the plot of the potential energy (which for a breather varies with time) and in fig 3 we plot the time 
dependence of the position of the local maximum of the energy density in the region of x < (which 
for a breather corresponds, alternatively to a kink and antikink). Of course, we could have given 
similar plot of the complementary kinks/antikinks which exist for x > 0. Our plots demonstrate that 
our initial configuration has turned into a breather, and as one can easily check - the frequency of the 
breather (which is related to its 'extent') is determined by the initial configuration (ie its d). 

We see that from t ~ 100 the field is essentially of the breather type. In fact - this can be checked 
by looking at the field itself at various times. Looking at the plots in fig. 3 we can evaluate the periods 
of oscillations of the breathers. We find that for d = 5 this period is approximately 436.29 giving 
to = 0.0134, for d = 3 its is T '^ 63.28 - giving uj = 0.094 and for d = 4 (not shown here) - it is 
T ~ 171.61 corresponding to w ~ 0.0366. As total energy is given by 



E = 2EoVl-u;^, (9) 

where Eq is the energy of one kink (in our case normalised to 4) we see that the corresponding values 
of energies are 

Ed=5 r^ 7.99927, Ed=4 ~ 7.9948, Ed=3 -- 7.9605 (10) 

which are in excellent agreement with our numerical results. 
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Figure 1: Total energies of field configurations seen for n = 2, started initially with a kink and an 
antikink separated hy a) d = 3.0, h) d = 5, c) d = 8. The little 'blips' are the numerical artifacts of 
calculation of energies when the breather is very thin. 
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Figure 2: Potential energy as function of time for n = 2 for the case of initial separation d = 5 
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Figure 3: Trajectory of a kink or antikink (in x < 0) for the case of n = 2 when the initial configuration 
had d = 5. 

Hence, the obtained field configurations are indeed, those of corresponding breathers. 
The same method can be used for the generalised models of Bazeia et al for other values of n. In 
the next section we look in detail at the case of n = 1. 



3 n = 1 model 

In this case the potential takes the form 



1 sm'^{2(f) 



fill 



2(1 + sin 99)2 

This potential, due to the factor (l + sinc/?)^ in the denominator, clearly, has a very different behaviour 
for 99 > and 99 < 0. Thus if we are to have a breather-like behaviour the field will vary over a different 
range of values for 99 > and 99 < 0. Of course, the initial kink and antikink can be chosen so that 
99 >0. 

We have performed several simulations (starting with different values of d). In fig 4. we present the 
results of the simulations corresponding to d = 1, 3, 5 and 9.5 (shown in fig 4a,b,c and d, respectively). 
We have found that all such configurations gradually lose energy but that this loss is very slow. 
Looking in detail at the plot of fig. 4 we note that this loss depends a lot on the initial separation (ie 
d). The energy drops slowly then almost settles and then drops again after which it gradually settles 
again etc. The field itself evolves very much like a breather (with unequal movement 'up' and 'down' 
- for the reasons mentioned above). However, the amplitude of this movement gradually decreases. 
This we can see from individual plots of the fields or from the time evolution of the position of the 
kink (or antikink) in x > 0. This position can be determined either by looking at the values of x at 
which the field takes a specific value, or even better, at which the energy density has a local maximum. 
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Figure 4: Total energies of the system for n = 1, started with a) d = 1, h) d = 3, c) d = 5 and d) 
d = 9.5. 



In fig 5. we present the plots of the time dependence of the position of the maximum of the energy 
density in x > 0. We note the above mentioned decrease of the amphtude, albeit a very slow one. In 
fact the decrease gets slower and slower suggesting that even though breathers probably do not exist 
the model possesses 'breather-like' configurations which are very long lived. The actual details depend 
a lot on the distance between the kink and the antikink of the initial configuration. 

Our results suggest that, perhaps, the nonexistence of breathers in this model is associated with 
the asymmetry of the potential as x ^^ —x. Hence, it is interesting to look in detail at the case of 
n = 3 (where this effect is reduced) and at n = 4 where the potential is symmetric. This we do in the 
next two sections. 
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Figure 5: Positions of the kink (or antikink) in x > for the simulations of n = 1 case in which the 
initial configurations were started with a) d = 1, h) d = 3, c) d = 5 and d) d = 9.5. 



4 n = 3 model 

First we consider n = 3. This time the potential takes the form 

Like in the previous case - the behaviour for ip > and c/? < is different; however, as compared 
to the n = 1 case the asymmetry of the potential is somewhat reduced. 

Again, we have performed many simulations for various values of d. In fig 6. we present plots of 
the energies for three such simulations (corresponding to d = 2,4 and 7). We note that, as before, 
the kinks and antikinks move towards each other and then interact with each other sending out a lot 
of energy (much higher percentage than in the n = 1 case discussed before). It seems that this time 
the fields annihilate quite quickly and do not form a long-living 'breather' - or if they do this breather 
is quite broad. To test this last idea we have decided to rerun one of the cases considered before , 
namely d = 4, with larger dx - so that our lattice is larger, namely —75 < x < 75; we show it as fig. 
6d. We do not see much difference - suggesting that the effects we have seen are genuine; ie that the 
field decays much more quickly than in the n = 1 case. 

5 n = 4 model 

Finally we have also looked at the n = 4 case. This time the potential is given by 

F((^) = ^sin2(2(^)(l+sinV)' (13) 

and, like the Sine Gordon case, the potential exibits symmetry x ^^ —x. 

However, performing the simulations in this case we have quickly discovered that this symmetry 
does not matter; all cases involved quick annihilation leading to pure radiation. In fig 7. we present 
the plots of the total energy corresponding to the cases of d = 6 and d = 10. We see a slow change 
of energy (when kink and antikink approach each other) followed by a rapid annihilation. We have 
looked at many other cases - but they produced qualitatively the same results (varying, of course, in 
the time required for the original approach). Hence, we feel that the n = 4 model is unlikely to possess 
breather like solutions. 

6 More general n 

We have also considered breathers for other values of n, In all the cases we have studied the kink 
and antikink approached each other and then annihilated straight away so it is clear that from all 
the models only the n = 2 possesses real breathers and for n = 1 the fields evolve into 'breather-like' 
configurations which do not appear to be stable but seem to be quasi-stable {ie are long-lived). 

The results mentioned above were obtained in the conventional discretisation (ie the discretisation 
in which the derivatives were represented by the symmetric differences of the fields etc) . The question 
then arises as to whether our results are 'robust' with respect to the change of discretisation. To check 
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Figure 6: Total energies of the system for n = 3, started with a) d = 2, b) d = 4, c) d = 7 and d) 
d = 4 but with larger dx 
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Figure 7: Total energy of the system started with a) d = 6, h) d = 10. 

that we can change the parameters of the discretisation and rerun our numerical simulations or, go 
further and consider completely different discretisations. Of these, perhaps the most interesting is the 
discretisation of Speight and Ward [3] as it preserves the Bogomolnyi bound [5]. Hence in the next 
section we discuss such discretisations applied to our case. 

7 Topological discretisation 



Consider the first class of Bazeia etc all models and look at their Lagrangian (1) and (2). Take the 
static part and rewrite (1) and (2) as 

2 



,^ ,2 /2A 2Asin"+V 

L = [Ox^j + — tan(/9 H 

\n n cos 99 



(14) 



2A 

Ox^ + — tan 99 
n 



2 A sin 



n+l 



95 



^ /2A 2Asin"+V' 

20x'^ — tanc^ + 

\ n n cos 99 



n cos Lp I 

However, when we calculate the total energy we integrate all terms above over x and we spot that the 
last term is a total divergence as we can write 

2Asin"+V^ 



^ / 2A 

20x'^ — tan 99 + 

\n n cos 99 



dxF, 



(15) 



where the explicit form of F can be easily calculated for each n. 

Hence after the integration over x the term (15) gives only the boundary contribution F{x = 
00) — F{x = -co). So, in the continuum limit, we have the Bogomolny equation for 99 which is given 

by 

2Asin"+V\ 



/2A 

Ux^ = — — tan 99 + 

\n n cos 99 



(16) 
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and the energy of all fields for which F{x = oo) — F{x = — oo) is given is smallest for the fields which 
satisfy (16). Thus the static kink is given by the solution of (16). As this is a first order equation for 
(f and is invariant under translations in x, its solution is unique, up to the constant of integration xq, 
which comes through ip = ip{x — xq) where xq is the position of the kink. 

The discretisation of Speight and Ward [3] is based on the observation that one can perform the 
discretisation of the Lagrangian (14) so that (15) still holds. Notice that this is possible as the two 
terms on the left hand side of (15) are the two terms that appear in the Langrangian (1) and (2). 
Hence with such a discretisation the minimum field is the field which solves the discrete Bogomolny' 
equation (again involving the same two terms) and its static solution is exact on the lattice. Thus 
such a field does not experience the Peirels Nabarro potential and, as shown by Speight and Ward 
in their paper [3] and in many subsequent papers [6], the numerical work with such fields has fewer 
numerical artifacts. 

So, here we follow the ideas of Speight and Ward and consider topological discretisation of (1) and 
(2). However, a little thought shows that this is not so simple. For n ^ 2 the expression for (2) is 
quite complicated and so the problem is not very simple. 

Notice that our expressions for F (15) are very different for n odd and n even. For n even we have 
(ignoring overall factors "* 



n' 



etc. 



sin Lo 1 

n = 2, V^tan^p ^, -^ F = -sin2(^, (17) 

cos 99 2 

Sin. LO t \ 

n = 4, F^tanw -, -^ F = -sin 209 sin4w, (18) 

cos ip 4 8 

sinV ^1 1 . 29 

n = 6, l/~tano9 , — > F= — sin 609 sin4o9 sin 209, (19) 

cos 99 32 4 32 

However for n odd we have 

sm (i9 

n = 1, y~tanci9 , -^ F = — log(l + sinw) + sinw, (20) 

cos 99 

■ 4 -1 

n = 3, F~tanv9- ^^^^, -^ F = - sin 4^9, (21) 

cos 99 8 32 

n = 5, V ^ tan 99 , — > F = — + sin Qip — — sin Atp — — sin 2^9, (22) 

cosc^ 32 192 64 64 

etc. 

Hence for n even we have (relatively) simple expressions for F involving sums of sin 2k'p> where k 
takes values up /c = n while for n odd we have in addition \og{(p). 

For n = 2 Speight and Ward took (for x = i dx) 

d,^ ^ [sin((c^(i + 1) _ ^(i _ 1))] _i_ (23) 

and 

sin(2v9) -^ sin(v9(i + 1) + ip{i - 1)) (24) 
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Both expressions have correct limits and the Bogomolny equation is satisfied as 

sin [ip{i + 1) + (p{i - 1)] sin [(p{i + 1) - (p{i - 1)] = - [cos(2^(i + 1)) - cos{2ip{i - 1))] . (25) 

Hence, when we integrate over x {ie sum over i), all the terms cancel except for the boundary terms. 
We can repeat this procedure for, say, n = 4 by taking again 

d^if ^ [sin((^(i + l)-c^(i_l))]_i_. (26) 

For the terms in the potential we can take 

sm{2ip) -^ sm{ip{i + 1) + (p{i - 1)) (27) 

and 

sin(499) -^ sm[ip{i + 1) + 2ip{i) + ip{i - 1)]. (28) 

Then the trigonometric identities mentioned above show that that this discretisation is again 
topological (we have the cancelation of terms involving cos(2c/?(i + 1)) and cos{2ip{i — 1)) as before and 
also of cos[2{ip(i + 1) + (p{i)] and cos[2{ip{i) + (p{i — 1)]). 

In the n = 2 case studied by Speight and Ward the Bogomolny equation (16) involves just two sin 
functions and is easy to solve. In our case (ie for n = 4) the equation involves three sin functions (of 
3 different arguments) and so is much harder to solve. This procedure is even harder to apply to the 
case of n = 6 where we have even more terms. And it seems almost impossible to generalise to the 
n = odd cases. And the lowest of these cases ie n = 1 and n = 3 are the most interesting from our 
point of view. 

Hence we have to look at the discretisation problem differently. 

7.1 A modified approach 

Let us look at the case of n = 1. In this case the potential is proportional to 

sin (fi , , 

F ~ tanv? -. (29) 

cos 99 

Naive integration, mentioned above, gives us 

F = -log(cos(^) - 2Arcth{tan^) + sin^?. (30) 

However, a few trigonometric identities tell us that F is also given by 

F = sin 99 — log(l + sin(^). (31) 

This is simpler but, at first sight, does not help us in our discretisation procedure. However, as we 
want to have the discrete terms in the Bogomolny bound cancel it is clear we would be happy if our 
discretisation gave us 

F -^ sm{ip{i + 1)) - sm{ip{i - 1)) - log(l + sm[ip{i + 1)]) + log(l + sm[ip{i - 1)]) (32) 

12 



This suggests to us that for our discrete potential we can take 

V{ip{i + 1), ^{i - 1)) = 2 dx . f -—- (33) 

It is easy to check that this expression has the correct continuum limit giving (29). 
Hence our discrete Bogomolny equation becomes 

K{[sin{{ip{i + 1) - ip{i - l))]f = sin((^(i+l))-sin(v9(i-l)) - log(l+sin[v9(i+l)])+log(l+sin[(^(i-l)]). 

(34) 

We did not succeed to solve this equation analytically and solved it numerically. Then we used 

it to construct the initial configuration which we then evolved using the discrete equations involving 

the expressions above. Unfortunately, they reproduced the results obtained before (using the usual 

discretisation). 

8 Final Comments and Further Work 

In a way our results are disappointing; the only time when we generated breather solutions corresponds 
to the Sine Gordon case, where these results are well known. In all other cases the initial configurations 
involving one kink and one anti-kink attracted, came towards each other and then annihilated. Of all 
our cases the n = 1 model was the nearest to possessing breathers; in its case the annihilation was 
quite slow and we have found that the model possesses 'breather like' long living states. In fact, such 
states have lived amazing long times; as is clear from fig 5, even after thousands of oscillations the 
configurations continue resembling a breather. 

Of course, the nonexistence of breathers in all models apart from n = 2 is likely to be related to 
the integrability of the Sine Gordon model, and non-integrability of the models for n ^ 2. On the 
other hand the fact that the n = 1 model has long-living states that resemble breathers suggests that 
this model is not far off from being integrable, ie is 'almost integrable'. We have no definition of this 
concept and we currently thinking of how to define it. 

We have also looked at two kink configurations for even n (at odd n such configurations would 
have infinite energy). In the Sine-Gordon models kinks repel (although some people argues that they 
go through each other). The preliminary results show that the same is true for other values of even n. 
In fig. 8 we present the trajectories obtained by us in two simulations (both started at d = 9 and with 
the same velocity v = 0.5. In fig 8a we present the trajectory for the Sine Gordon case n = 2 and in fig 
8b - for n = 4. We see very little difference although the repulsion is stronger in the n = 4 case. We 
observed the same results for other values of v and for different values of d. In fig 9 we show similar 
results for the case of !> = (this time started at a smaller d - at d = 9 but the repulsion is much 
weaker). Clearly we see repulsion in both cases. In our simulations we took the initial conditions as 
the superpositions of two kinks with velocity ±v in their centre of mass. 

Of course, in the Sine Gordon model one has analytical expressions for moving kinks; and it is not 
clear who to relate the numerical results to the analytical ones. 

We plan to look at this in more detail in our further work. 
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Figure 8: Trajectories of two kinks started with d = 9 and with velocities v = 0.5 in two models of 
even n: a) n = 2, b) n = 4. 
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Figure 9: Trajectories of two kinks started with d = 6.5 and at rest (ie with f = 0) in two models of 
even n: a) n = 2, b) n = 4. 
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